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In this work we study the properties of a two-subband quasi-two-dimensional electron system 
, in a strong magnetic field when the electron filling factor is equal to four. When the cyclotron 

' energy is close to the intersubband splitting the system can be mapped onto a four-level electron 

I system with an effective filling factor of two. The ground state is either a ferromagnetic state 

or a spin-singlet state, depending on the values of the inter-level splitting and Zeeman energy. 
Q ' The boundaries between these phases are strongly influenced by the inter-electron interaction. A 

, significant exchange-mediated enhancement of the excitation gap results in the suppression of the 

■ electron-phonon interaction. The rate of absorption of non-equilibrium phonons is calculated as a 

' function of Zeeman energy and inter-subband splitting. The phonon absorption rate has two peaks 

. as a function of intersubband splitting and has a step-like structure as a function of Zeeman energy. 
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During the last two decades the behavior of a two-dimensional (2D) electron system in a quantizing magnetic field 
attracted a lot of theoretical and experimental interest. The peculiar properties of this system result from its effective 



^ . zero dimensionality, because the magnetic field normal to the 2D layer quenches the electron kinetic energy to a 
constant. As a result the Hamiltonian of the system consists of the interaction part only, and the system undergoes 
the transition into the incompressible state at certain magnetic field values. The formation of the incompressible 
■ states results in the Fractional Quantum Hall Effect - one of the most interesting phenomena discovered in the 

a ' strongly correlated electron systcms.El^cl The only typical energy of the system is the characteristic electron-electron 
Coulomb interaction energy, which is of the order of ec — / kI, where I is the magnetic length and k, is the dielectric 
constant. Introduction of a new degree of freedom into this system results into new interesting phenomena when the 
, characteristic energy of the new degree of freedom becomes smaller than the Coulomb energy, ec- The electron spin is 
Q one of the possible degrees of freedom. If the Zeeman energy is small enough then the trajisition from spin-polarized 
O \ to spin-unpolarized ground state occurs at electron filling factor vq = 2/m, where m is odd.u At filling factor i/q = 1/m 
a new type of charged topological excitations, skyrmions, appears. 
' One can also introduce the degree of freedom in the direction normal to the 2D layer. This can be realized in a 
^ ■ double-layer system. In this case the phase diagram is driven by the interplay between the electron-electron Coulomb 
OO ■ energy, the tunnelling energy between two layers and the Zeeman energy. It was found that at i^o = 2 the double-layer 
system can be found inlhree different phases: fully spin-polarized ferromagnetic state, spin-singlet state and canted 
' antiferromagnetic state.El 

^ , In this paper we consider the system which is similar to a double-layer system but has a different geometry, which 
I ■ provides more experimental possibilities to change the parameters of the system. Namely, we study a single hetero- 
junction in which the second degree of freedom is introduced by the size quantization in the direction perpendicular 
"Y^ to the 2D layer. If the electron filling factor is greater than two and the energy A = A12 — hujc is smaller than 
d ' the Coulomb energy, then the second subband has to be taken into account. Here huJc is the cyclotron energy and 
^ . A12 is the inter-subband splitting. In this case the second Landau level (Landau level number n = 1) of the first 
I subband is close in energy to the first Landau level (n = 0) of the second subband. Since the intersubband splitting 
A12 has a weak dependence on the magnetic field, the separation A can be changed by changing magnetic field B. 
The influence of the second subband on the optical properties of the magnetically-quantized quasi-2D systerns was 
observed in the different types of experiments, e.g., in extrinsic radiative recombination magnetospectroscopyB and 
in the optically-detected cyclotron resonance in tilted magnetic field.l3 It was proposed in Ref. ^ that the intersection 
between the levels can be also observed in the magnetic field dependence of the phonon-mediated conductivity of the 
. system. At filling factor slightly greater than two, the dissipative conductivity should reveal a double-peak structure 
$H as a functiorLX)|Lmagnetic filed. A similar double-peak structure should be observed in the other phonon spectroscopy 
5^ , experiments.LrEd With increasing electron density the repulsion between the levels due to electron-electron exchange 
interaction opeiia^a gap in the inter-level excitation spectra. As a result the double-peak structure transforms into a 
single-peak oncJlj 

In what follows we study the properties of the two-subband system when the electron filling factor is equal to four, 
= 4. We assume that the Coulomb energy is much smaller than the cyclotron energy. In this case the completely 
occupied lowest Landau level of the first subband can be considered as a non-dynamical background. We will be 
interested only in a subsystem consisting of the second Landau level (n = 1) of the first subband and the first Landau 
level {n = 0) of the second subband, the spin degeneracy of both levels is lifted by the Zeeman splitting A^. This 
four-level subsystem is schematically shown in Fig. |l|. For vq = A the effective filling factor v of this subsystem is 



1 



equal to two, v = 2. Under this mapping the system becomes similar to a double-layer system with the total filling 
factor V = 2 and for its analysis we follow the method of Ref. ||. 

For non-interacting electrons we can distinguish three cases: 1) if A > the states '1' and '2' are occupied, and 
the ground state is a spin singlet state; 2) if — A^ < A < A^ the states '1' and '3' are occupied, and the ground 
state is a ferromagnetic state; 3) if A < — A^ the states '3' and '4' are occupied, and the ground state is again a spin 
singlet state. The boundaries between these phases are given by the equations: A = A^ and A — — A2. When the 
electron-electron interaction is taken into account the transitions between spin singlet and ferromagnetic states can 
occur through a new phase, for example through a canted antiferromagnetic phase.B Our calculations show that in 
the present system there is no intermediate phase and the transitions from the singlet to the ferromagnetic state are 
sharp, like in the non-interacting case. 

In what follows we use the Coulomb energy ec as the unit of energy and magnetic length / as the unit of length. 
In the Landau gauge with the vector potential A = (0, Bx, 0) the eigenstates of a single-electron Hamiltonian are 
characterized by the Landau level number, n; y-component of the momentum, fcj,; the electron subband number, 
/i = 1 or /i = 2, and the z-projection of the electron spin, a = 2Sz = ±1: 

^iky y 

■lpn,ky,t,,aix,y,z) ^ Xliiz)^a^=(t>nix - ky) , (1) 

where Xfii^) is the envelope wavefunction of the fith subband; ^o- is the spin part of the wavefunction; (j)nix) is the 
nth harmonic oscillator function. 

The Hamiltonian of the interacting electron system is 

2" 2^ ^ '•^)^'nkytia'~''ril<:yt^a 
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where and Cn,k,fia are the creation and annihilation operators of the electron in the state ipn,ky,fi,a- In Eq- (§) 

we use the notations: 

= ^^..u.3.2(9)G„,„,(9)G„3n2(r)exp (^-^) , (3) 

whereEl 

Q = Qx + i<ly, Q = \q\ and is a generalized Laguerre polynomial. The modification of the Coulomb interaction due 
to the finite extent of the electron wavefunctions in z-direction is given by: 

POO POO 

Fnu^2^3{<l) ^ j / dzidz2e^'^'^''^~''^'^Xii{zi)XiAzi)Xi2{z2)Xi3{z2) ■ 
Jo Jo 



We use the Fang-Howard approximationlli for the envelope wavefunctions of the electrons in the first and the second 
subbands: 

xi{zj = Y y^'^^p V 2 

X2W-/f^(^-^)exp(^-i6z 

We follow the standard Hartree-Fock method, assuming a non-zero average of (^C^^ fi^cri^k2,p.2<y2'j o'^^r the ground 
state. The specific feature of the given problem, which differs it from the double-layer system, is that the non-zero 
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paring in (||) can occur between the states with the different value of ky, because in our case the two-level system is 
formed by the different Landau levels.cJ For each value oi ky = k we introduce the new wavef unctions, which are the 
eigenf unctions of the Hartree-Fock Hamiltonian. The creation and annihilation operators corresponding to the new 
wavefunctions are a^^, Ofe^i, where i = 1, 2, 3 or 4. These functions are related to the original ones by the matrix ^i^y'. 

L'fe+Q,2,-l — Z^i=i liiO'k,! 
Ck+Q.2,1 = '^i=lliiO-k,i ■ 



The average of the introduced functions over the ground state is {a'^iak,jj = Sij^Sn + Si2)- This means that only the 
states with the lowest energies (i = 1 and i = 2) are occupied. Substituting Eq.(^ into the Hamiltonian (|^) we obtain 
the Hartree-Fock Hamiltonian in the form of a four-by-four matrix. In the basis {Ck,i,-i, Ck,i,i, Ck+Q,2,-i, C'fc+Q,2,i) 
the elements of this matrix are 

H^/ = -^(A + A,) - /3ne}l - P33e% 
Hg^ = -^(A - A,) - /322e}} - 



= -/3i3(<S(o)+<-;(Q); 

HF _ T/(0101).^^ _ a .t/(0011) 



ttHF 



where /3y = 7ii7ji + 7i27j2, 

^uS^^ ^ dkPnnnjAk) \G^,^M\ \G^,^^k)\ J|,,+,3_,,_,,| (fcg) exp (^-^^ , (5) 



ili2 „ t/(h*2 il 22) /r)\ 

Jm is the Bessel function of the mth order, and e*^^^^ is the exchange energy of an electron in the iith Landau level 
of the jith subband interacting with electrons of the same spin in the filled i2th Landau level of the j2th. subband. 

Taking into account that the eigenvectors of the Hartree-Fock matrix ak^i are related to Ck,^,cr by the matrix 7^ 
(Eq. (^), we obtain the self-consistent system of equations for 7^ . The parameter Q can be found by minimizing the 
Hartree-Fock energy. We have found the solution of the system of equations for 7jj numerically. For all values of A 
and Az the ground state of the system is one of the non- interacting phases (spin singlet or ferromagnetic) , described 
above. The transitions between the phases are sharp. The interaction modifies the phase boundaries only: 

A = A^ + \{4l-e\\)+eZ-ef,+ell = A. + eH = AH , (6) 
A = -A,-fi(e^^-.}})+^?2-^?}-^?^ = -A. + 6,^A^ , (7) 

where A// and A^ denote the higher and lower values of splitting A at the boundary of the ferromagnetic phase, 
see Fig. ||(a). In Fig. ||(a) the phase diagram is shown for 6 = 1. The phases 5*1 and S2 are the spin-singlet phases 
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and the phase F is the ferromagnetic phase. The single-electron states '1' & '2', '1' & '3', and '3' & '4' (Fig. [|) are 
occupied in the phases 5*1, F, and 52, respectively. The electron-electron interaction results in the splitting of the 
boundaries between different phases, even when the Zeeman energy is equal to zero. For = the splitting is equal 
to €h — — '^^127 where en and cl are introduced in Eqs. and (0). In Fig. ^(b) the values of the lower (ei,) 
and upper (en) phase boundaries at A2 — are shown as the function of the Fang- Howard parameter b, where 1/6 
is proportional to the characteristic width of the heterojunction in units of magnetic length. Both values en and 
decrease with increasing b. 

We are mostly interested in the case when the intersubband splitting is close to the cyclotron energy. This gives the 
restriction of the values of the parameter b. This parameter should be close to the inverse magnetic length. However, 
the exact relation between b and / for A = depends on the actual value of b. Therefore we found it appropriate to 
calculate the 6-dependencies of the main parameters of the system. 

The two-subband system can be studied experimentally by changing magnetic field when keeping the fixed value of 
the filling factor. If the magnetic field B is normal to the heterojunction, then both —A and A^ increase linearly with 
increasing B. In this case only the transition from the phase to the phase F can be observed. Another possibility 
is to make the experiments in the tilted magnetic field with the fixed value of its normal component. In this case A 
is a constant, because it is proportional to the normal component of the magnetic field, but A^ changes linearly with 
the total magnetic field. Depending on the initial value of A, one should observe the transition from 5*1 to F or from 
^2 to F phases. However, there is a finite region of A in which no transitions can be found and the system is always 
in the phase F. This region has a non-zero range even for zero A^ when its width is 2el2- 1-1 1—1 

Now we shall discuss the possibility to detect the different phases by the acoustic-phonon spectroscopy.Eitil We 
consider the absorption of non-equilibrium phonons_by the discussed quasi-2D system. The electron-phonon interaction 
is proportional to the electron density opcrator.ll3 During the act of the phonon absorption the electron system 
undergoes the transition from the ground state to the excited state, corresponding to the density fluctuations. When 
the system is in the ferromagnetic phase (F), there are spin reversal excitations only. In this case the phonon 
absorption is forbidden. The collective exoitation spectra in the singlet phases and ^2 can be found from the poles 
of the density-density correlation functionEj, which are given by the expressions: 

E{q) = A + 2F°2\°2^ (q) - V??^ (q) - e',l - e<il + + e?} (8) 

for the phase 6*1 and 

E{q) = - A + 2F°2\°2^ (q) - VZ\' (<Z) - e',l - 41 + ell + 41 (9) 

for the phase ^2, where Eqs. (||) and (||) were used. One can rewrite Eqs. (^) and in the form, which can describe 
both cases simultaneously: 

Eiq) = \SA\ + A, + 2V,%\%\q) - V,%\\q) + \ {e^ + ell) , (10) 

where SA — A ~ Ah for the phase Si and ^A = —A + Al for the phase ^2. 

The energy spectra E{q) are shown in Fig. ^(a) for SA = and A^ = 0, for several different values of parameter 
6. One can see that there is a finite gap for all values of momentum q. This means that the gap exists for any 
value of A and A^. With increasing the spreading of the electron wavefunction in z direction (decreasing b) the gap 
becomes smaller, which results from the decreasing of the effective 2D electron-electron interaction. The existence 
of the interaction-induced finite excitation gap in the two-level system for any value of the inter-level splitting was 
demonstrated for the system with filling factor v — 1 in Ref. The effect of the strong renormalization of the 
excitation energy due to the iiiteij-«lectron interaction is also known for spin excitations as an interaction enhancement 



(11) 



of the g-factor in 2D systems Jla~l2i 

The electron-phonon interaction Hamiltonian has the form: 



3,Q 



where the isotropic Debye approximation is used with the linear dependence of the phonon frequency on the wave 
vector: ujj{K) = SjK, sj being speed of sound; j is labeling the phonon modes, j = 1 for longitudinal mode and 
j = 2,3 for two transverse modes. In Eq. (|l^) we labeled the three dimensional (3D) phonon wave vector by the 
capital letter, Q, and its projections by the corresponding small letters, Q = {q,qz)- The creation operator of a 
phonon in the jth mode is denoted as d^, V is a normalization volume, p{q) is the two-dimensional electron density 
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operator, Mj{Q) is the matrix elemfiat of the electron-phonon interaction, which is determined by the deformation 
potential and piezoelectric coupling:tJ 



2/90 



rjQxQy^j.z + QyQ z^j,x + Q zQx£,j.y , ^ /^^ 

-P g2 ' ^■^OUi • Q) 



(12) 



where po is the mass density, (3 and Sq are the parameters of piezoelectric and deformation potential couplings, is 
the polarization vector of the jth phonon mode. The parameters of GaAs arc used in our calculations. 

Because the phonon-assisted transitions are only allowed between the different subbands, the form factor Z{qz) is 

Z{q,) = j dze'''''xi{^)X2{z) . (13) 

At low teBwerature, fc^T <C E{q), the rate of absorption of non-equilibrium phonons can be found from the 
expression :E5tZl 

_2tt f dQ 



^SiEiq)-s,Q)n,{Q)\M,{Q)Z{qz)\\oi{q) , (W) 



where Roi{q) = (g'^/2) exp(— g^/2), nj{Q) is the phonon distribution function. 

One can see how the appearance of the excitation gap suppresses the ph onon ab sorption. The phonon absorption 
process requires the conservation of energy and in-plane momentum, Sj y^q^ + q-^ = E{q), where the characteristic 
excitation gap, E{q), is greater than 0.1 for 6^1 (see Fig. ^a)). The GaAs speed of sound in units of lec/^ is about 
0.03 for longitudinal and 0.02 for transverse phonons. Therefore, the large value of \/ q'^ + q^ is required to satisfy the 
energy conservation law. However, the factors i?oi('z) and Z{qz) in Eq. ( p^ ) make the phonon absorption to be small 
at g > 1 or > 6. Thus, the phonon absorption becomes strongly suppressed comparing to the non-interacting case, 
in which the excitation energy A is allowed to have any value. For the non-interacting electrons the absorption rate 
has two maxima at |A| ^ hsj/lB 

From Eq. (|l^ ) we calculate the normalized phonon absorption rate, ujq ~ ujabs/n{E{qo)), where go is a characteristic 
momentum, here we use q^ = 1. In Fig. |^(b) the rate is shown as a function of A for = and different values 
of h. The absorption rate has a two-peak structure. The gap between the peaks corresponds to the phase F of the 
system, in which the phonons can not be absorbed. The absorption rates for 6 = 1 and 6 = 1.5 are multiplied by 
10. With increasing h the collective excitation gap increases (see Fig. ||(a)), which results in a rapid decrease of the 
absorption rate. One can see this tendency for h = 0.5 and 6=1. However, with increasing b the phonons with 
the larger values of qz can be absorbed, which tends to increase the absorption rate. The latter effect becomes more 
pronounced for large h when the excitation gap has a weaker dependence on h (see Fig. |^(a)). The competition 
between these two effects results in the slight increase of the absorption rate for b — 1.5. 

In the tilted magnetic field experiments, when the Zeeman energy A2 is changed while the inter- level splitting A 
is kept fixed, the absorption rate remains constant. For the singlet phase 5*1 or ^2 the absorption rate is given by 
Eqs. (H), (^), ( p^ and for the ferromagnetic phase F the absorption rate is zero. The dependence of ujq on A2 is 
shown schematically in Fig. |^(c), where the critical value of Zeeman splitting. A*, corresponds to Si-F or S2-F phase 
transition. 

In conclusion, we considered the two-subband electron system with the total filling factor four. The system was 
mapped onto the four-level electron system with the effective filling factor two. The separation between the levels 
belonging to the different subbands is proportional to the cyclotron energy and can be changed by changing the 
magnetic field. The system exists in one of the three different phases: one is ferromagnetic and two are spin-singlet 
phases. The electron-electron interaction does not create new phases. However, it renormalizes the phase boundaries. 
The excitation spectrum has the gap for any value of the inter- level splitting A and Zeeman splitting A2 . This results 
in the strong suppression of the electron-phonon interaction. The rate of the phonon absorption by the considered 
quasi-2D electron system has a double-peak structure as a function of level splitting and a step-like structure as a 
function of Zeeman splitting. 

The authors acknowledge the support by the UK EPSRC. 



* E-mail address: m.e.portnoi@ex.ac.uk 



5 



^ R.E. Prange and S.M. Girvin, The Quantum Hall Effect (Springier, Heidelberg, 1987). 

^ T. Chakraborty and P. Pietilainen, The Fractional Quantum Hall Effect (Springer, New York, 1995), 2nd ed. 
^ Perspectives in Quantum Hall Effects, edited by S. Das Sarma and A. Pinczuk (Wiley, New York, 1997). 

T. Chakraborty, Adv. Phys. 49, 959 (2000). 
^ S. Das Sarma, S. Sachdev, and L. Zheng, Phys. Rev. B 58, 4672 (1998). 
® A.S. Plant, I.V. Kukushkin, K.v. Klitzing, and K. Ploog Phys. Rev. B 42, 5744 (1990). 

D.D. Arnone, T.P. Marlow, C.L. Foden, E. H. Linfield, D. A. Ritchie, and M. Pepper, Phys. Rev. B 56, R4340 (1997). 
* V.M. Apalkov and M.E. Portnoi, Physica E, to be published; V.N. Golovach and M.E. Portnoi, unpublished. 
® L.J. Challis and A.J. Kent, in Aspects of Semiconductor Nanostructures, edited by G. Bauer, F. Kuchar and H. Heinrich 

(Springer, Heidelberg, 1992). 

^° C.J. Mellor, R.H. Eyles, J.E. Digby, A.J. Kent, K.A. Benedict, L.J. Challis, M. Henini, C.T. Foxon, and J.J. Harris, Phys. 
Rev. Lett. 74, 2339 (1995). 

" U. Zeitler, A.M. Devitt, J.E. Digby, C.J. Mellor. A.J. Kent, K .A. Benedict, and T. Cheng, Physica B 249-251, 49 (1998). 

V.M. Apalkov and M.E. Portnoi, preprint cond-mat/0111377 
" C. Kallin and B.I. Halperin, Phys. Rev. B 30, 5655 (1984). 

T. Ando, A.B. Fowler, and F. Stern, Rev. Mod. Phys. 54 , 437 (1982). 
" G.F. Giuliani and J.J Quinn, Phys. Rev. B 31, 6228 (1985). 

V.F. Gantmakher and Y.B. Levinson, Carrier Scattering in Metals and Semiconductors (North-Holland, Amsterdam, 1987). 
" K.A. Benedict, R.K. Hills, and C.J. Mellor, Phys. Rev. B 60, 10984 (1999). 

R.J. Nicholas, R.J. Hang, K.v. Klitzing, and G. Weimann, Phys. Rev. B 37, 1294 (1988). 
" B.B. Goldberg, D. Heiman, and A. Pinczuk, Phys. Rev. Lett. 63, 1102 (1989). 
^° A. Usher, R.J. Nicholas, J.J. Harris, and C.T. Foxon, Phys. Rev. B 41, 1129 (1990). 

A.P. Smith, A.H. MacDonald, and G. Gumbs, Phys. Rev. B 45, 8829 (1992). 



FIG. 1. The single-electron energy levels of the two-subband system. The levels '1' and '2' belong to the second Landau 
level {n = 1) of the first subband with spins 5*^ = 1/2 and Sz ~ —1/2, respectively. The levels '3' and '4' belong to the first 
Landau level (n = 0) of the second subband with spins Sz = 1/2 and Sz = —1/2, respectively. The splitting A is the interlevel 
splitting, Az is the Zeeman energy. 
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FIG. 2. (a) The phase diagram of the four-level system for 6=1. The phases Si and S2 are spin singlet phases, the phase 
f is a ferromagnetic phase, (b) The parameters e^r and (Eqs.(6)-(7)) of the phase diagram are plotted as the functions of 
b. All energies are in units of sc, b is in units of 1/Z. 




FIG. 3. (a) The excitation dispersion E{q) (Eq. (10)) for 5 A = = A^. The numbers near the lines show the values of the 
parameter b, where b and q are in units of 1/Z, E{q) is in units of ec- (b) The phonon absorption rate wo as a function of A for 
Az — 0. Solid, dotted and dashed lines are for b = 0.5, 1.0 and 1.5, respectively. The data for b — 1.0 and 1.5 are multiplied by 
10. The absorption rate uiq is in units of 10^'^ s~^, A is in units of ec- (c) Phonon absorption rate oJo as a function of Az for 
constant A. Here u!o and A^ are in arbitrary units. The critical Zeeman energy A* corresponds to the Si—fF or phase 
transitions. 
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